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Abstract: We investigate the phase transition in fuh QCD with two flavors of staggered 
fermions in presence of a constant abehan chromomagnetic field. We find that the critical 
temperature depends on the strength of the chromomagnetic field and that the deconfined 
phase extends to very low temperatures for strong enough fields. As in the case of zero 
external field, a single transition is detected, within statistical uncertainties, where both 
deconfinement and chiral symmetry restoration take place. We also find that the chiral 
condensate increases with the strength of the chromomagnetic field. 
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1. Introduction 

In previous studies [1,2] on the vacuum dynamics of pure non-abelian gauge theories we 
found that the deconfinement temperature depends on the strength of an external abelian 
chromomagnetic field. In particular we ascertained that the deconfinement temperature 
decreases when the strength of the applied field is increased and eventually goes to zero. 
It is not difficult to see here an analogy with the reversible Meissner effect in the case 
of ordinary superconductors (where the system goes to normal even at zero temperature 
if the magnetic field is strong enough) and therefore we referred to it as "vacuum color 
Meissner effect". We have also verified that the same effect is not present in the case of 
abelian gauge theories, so that it seems to be directly linked to the non-abelian nature of 
the gauge group. 

The dependence of the deconfinement temperature on applied external fields is surely 
linked to the dynamics underlying color confinement, therefore in our opinion, apart 
from possible phenomenological implications, such an effect could shed light on confine- 
ment /deconfinement mechanisms. 

On these basis we believe that it is important to test if the effect continues to hold 
and how it qualitatively changes when switching on fermionic degrees of freedom. One 
important aim of the present work is therefore to investigate the dependence of the decon- 
finement temperature on the strength of an external abelian chromomagnetic field in the 
case of full QCD with two flavors. 

A second important and relevant issue regards the relation between deconfinement 
and chiral symmetry restoration. As it is well known, the two phenomena appear to be 
coincident in ordinary QCD, while they are not so in different theories (like QCD with 
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adjoint fermions [3-5]). A simple explanation of this fact is not yet known and may be 
strictly linked to the very dynamics of color confinement. An important contribution 
towards a clear understanding of this phenomenon could be to study whether it is stable 
against the variation of external parameters: both theoretical and numerical studies [6-9] 
have been performed in that sense for the case of QCD in presence of a finite density 
of baryonic matter. In the present work we investigate the same issue for the case of an 
external field, i.e. we will ascertain whether deconfinement and chiral symmetry restoration 
do coincide also in presence of a constant chromomagnetic field. 

The paper is organized as follows. In Section 2 we review our method to treat a 
background field on the lattice. In Sections 3 and 4 we discuss numerical results and 
finally, in Section 5, we present our conclusions. 

2. External fields on the lattice 

In this section we review our method to study the dynamics of lattice gauge theories in pres- 
ence of background fields. In particular, we focus on the case of constant chromomagnetic 
fields. 

2.1 The method 

In Refs. [10, 11] we introduced a lattice gauge invariant effective action r[A'^^*] for an 
external background field A^^^: 

where Lt is the lattice size in time direction and A'^^^{x) is the continuum gauge potential 
of the external static background field. ^^[A^^*] is the lattice functional integral 

Z[A^^*] = [ VU e"^^ , (2.2) 

with Sw the standard pure gauge Wilson action. The functional integration is performed 
over the lattice links, but constraining the spatial links belonging to a given time slice (say 
xt = 0) to be 

Uk{x,xt = 0) = Ul^\x) , {k = 1,2,3) , (2.3) 

C/^^*(x) being the elementary parallel transports corresponding to the external continuum 
gauge potential A^^^{x) = A'^^{x)\a/2. Note that the temporal links are not constrained. 
Z[0] is defined analogously, but adopting a zero external field, i.e. with U^^{x) fixed to 
the identity element of the gauge group. 

In the case of a static background field which does not vanish at infinity we must also 
impose that, for each time slice xt 7^ 0, spatial links exiting from sites belonging to the 
spatial boundaries are fixed according to eq. (^.3|). In the continuum this last condition 
amounts to the requirement that fluctuations over the background field vanish at infinity. 
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The partition function defined in eq. ( |2.2D is also known as lattice Schrodinger func- 
tional [12, 13] and in the continuum corresponds to the Feynman kernel [14]. Note that, 
at variance with the usual formulation of the lattice Schrodinger functional [12,13], where 
a cylindrical geometry is adopted, our lattice has an hypertoroidal geometry, i.e. the first 
and the last time slice are identified and periodic boundary conditions are assumed in the 



time direction, so that the constraint given in eq. (22) should actually read Uk{x,Lt) = 
Uk{x,0) = Ul^^{x). With this prescription, Sw in eq. ( ^^ ) is allowed to be the standard 
Wilson action. 



The lattice effective action r[A^^*] defined by eq. (2T) is given in terms of the lattice 
Schrodinger functional, which is invariant for time-independent gauge transformation of 
the background field [12, 13], therefore it is gauge invariant too. It corresponds to the 
vacuum energy, So[^^^*]) in presence of the background field, measured with respect to the 
vacuum energy, -EqIO], with ^4'^^* = 

Pj^cxt] — ^ So[l^^*] -So[0]. (2.4) 

The relation above is true by letting the temporal lattice size ^ oo; on finite lattices 
this amounts to have Lf sufficiently large to single out the ground state contribution to the 
energy. 

For finite values of Lt, however, having adopted the prescription of periodic boundary 
conditions in time direction, the functional integral in eq. ( |2.2[) can be naturally interpreted 
as the thermal partition function Zq-[A^^^] [15] in presence of the background field A^^^, 
with the temperature given by T = l/(aLt). 

In this case the gauge invariant effective action in eq. (|2.1| ) is replaced by the the free 
energy functional defined as 



When the physical temperature is sent to zero the free energy functional reduces to the 
vacuum energy functional, eq. (l2.l|) . 

Let us now consider the extension of the above formalism to full QCD, i.e. including 
dynamical fermions, which is relevant for the present work. In presence of dynamical 
fermions the thermal partition functional becomes [16] 



Zt 



J JUk{Lt,x)=Uk(0,x)=U^^^(x) 

= / VUe-^"^ detM, (2.6) 

JUk{Lux)=Uk(0,x)=Ul^\x) 

where Sp is the fermion action and M is the fermionic matrix. The spatial links are still 
constrained to values corresponding to the external background field, whereas the fermionic 
fields are not constrained. The relevant quantity is still the free energy functional defined 



as in eq. (2.5) 



Actually, a direct numerical evaluation of the ratio of partition functions appearing in 
eq. ( |2.5| ) turns out to be quite difficult. Even if techniques have been developed recently 
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to deal with similar problems [17, 18], we adopt the more conventional strategy [10, 19] of 
computing instead a susceptibility of the free energy functional, in particular its derivative 
F' with respect to the inverse gauge coupling /?, which can be easily evaluated and is also 
more appropriate for the aim of the present study. F' is defined as 



F'(/3) 



dp 



J2 ^ReTrC/^,(x) 

\x,fj,<u 



(2.7) 



where the subscripts on the averages indicate the value of the external field. Only un- 



constrained plaquette are taken into account in the sum in eq. (2/7). Observing that 
^^^extj _ Q ^ _ j^g^y cventuahy obtain F[^'^''*] from F'lA'^^^] by numerical inte- 

gration: 



F[A 



exti 



F'[A<=^*](i/3' 



(2. 



2.2 A constant chromomagnetic field on the lattice 



Let us now define a static constant abelian chromomagnetic field on the lattice. In the 
continuum the gauge potential giving rise to a static constant abelian chromomagnetic 
field directed along spatial direction 3 and direction a in the color space can be written in 
the following form: 



At\x) 



{2S 



In SU(3) lattice gauge theory the constrained lattice links (see eq. (|2.3f) ) corresponding to 
the continuum gauge potential eq. (|2.9D are (choosing a = 3, i.e. abelian chromomagnetic 
field along direction 3 in color space) 



Uf^'^x) = t/f*(f) = 1, 

• agHx\ 



T\x) 



exp(i- 2 





exp( 











1 



(2.10) 



We will refer to this case as T3 abelian chromomagnetic field, which will be our choice in 
the present work. Of course it is possible to choose various alternatives, like an abelian 
field along the direction of the Tg generator, or along different combinations of and Tg. 

Since our lattice has the topology of a torus, the magnetic field turns out to be quan- 
tized 

In the following next will be used to parameterize the external field strength. 



2tt 

T ^ext 



next integer. 



(2.11) 



3. Numerical simulations and results 



We have studied full QCD dynamics with two flavors of staggered fermions in presence of 
a constant chromomagnetic field. The simulations have been performed on lattices 32^ x 8 
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and 64 x 32^ x 8. We used a slight modification of tlie standard HMC R-algorithm [20] 
for two degenerate flavors of staggered fermions with quark mass aniq = 0.075. According 
to our previous discussion, the links which are frozen are not evolved during the molecular 
dynamics trajectory and the corresponding conjugate momenta are set to zero. We have 
collected about 2000 thermalized trajectories for each value of /?. Each trajectory consists of 
125 molecular dynamics steps and has total length 1. The computer simulations have been 
performed using computer facilities at the INFN apeNEXT computing center in Rome. 

3.1 The critical coupling 

As is well known, the pure SU(3) gauge system undergoes a deconfinement phase transition 
at a given critical temperature and this happens even in the unquenched case (see Ref. [21] 
for an up-to-date review). In our earlier studies [1,2] we found that the critical coupling in 
pure non-abelian gauge theories is shifted towards lower values by immersing the system in 
a constant chromomagnetic background field: that means lower temperatures on lattices 
where the temporal extent in lattice units is kept constant (T = l/(Lta)). The main 
purpose of the present study is to verify if this effect survives in presence of dynamical 
fermions. We refer in the following to the constant abelian background field defined in 
Eqs. (^ and ( plO| ). 



In order to evaluate the critical gauge coupling we measure F'[^''^*] (eq. (^)), the 
derivative of the free energy with respect to the gauge coupling /3, as a function of /?. We 
found that F'[A''^*] displays a peak in the critical region where it can be parameterized as 

a2{Lt)[l3 - (3*{Lt)\^ + 1 

In figure 1 we show an example of F' measured for next = 1 on a 32^ x 8 lattice. In the 
same figure we display also the chiral condensate 

= {1- ^ Tr (3.2) 

and the numerical data point out that the peak in the derivative of the free energy corre- 
sponds to the drop in the chiral condensate, the latter is a signal of the transition leading 
to chiral symmetry restoration. 

In figure |2| we compare the derivative of the free energy with the absolute value of the 
Polyalov loop 

^=^E^Trn[/4(x4,f). (3.3) 

X X4 = l 

We can see that also in this case the rise of the Polyakov loop (expected at the deconfining 
phase transition) corresponds to the peak of the derivative of the free energy. Moreover, 
in figure ^ the derivative of the free energy is displayed together with the plaquette sus- 
ceptibility (susceptibility of the gauge action). It is evident from this figure that, within 
our statistical uncertainties, the peaks of the two quantities coincide. Similar results are 
obtained by looking at the susceptibilities of the Polyakov loop and of the chiral condensate. 
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Figure 1: The derivative of the free energy eq. (|2.7|) with respect to the gauge coupling (left axis 



blue circles), and the chiral condensate eq. (3.2) (right axis, red squares) versus (3. The vertical line 
represents the position of the peak in the derivative of the free energy. 
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Figure 2: The derivative of the free energy as in figure together with the absolute value of the 
Polyakov loop eq. (|3.3|). Vertical dotted line as in figure lu 



From the above arguments we may draw some partial conclusions: the critical coupling 
of the phase transition can be located by looking at the peak of the derivative of the free 
energy; moreover, as in the case of zero external field and within statistical uncertainties, 
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Figure 3: The derivative of the free energy as in figure |l| together with the susceptibility of the 
gauge action. 



a single transition seems to be present where both deconfinement and chiral symmetry 
restoration take place. 

It is worth to note that since the measurement of the derivative of the free energy is 
simply related to the measurement of the gauge plaquette we may have a good evaluation 
of the critical coupling with a relatively small sample of measurements. 

As told before our aim is to find if the critical coupling depends on the strength of the 
applied constant chromomagnetic field. To this purpose we have varied the strength of the 
external field by tuning up the parameter riext and we have searched for the phase transition 
signalled by the peak of the derivative of the free energy. We have found that indeed the 
critical coupling shifts towards lower values by increasing the external field strength. In 
figure 1^ we display the derivative of the free energy in correspondence of three values of 
next obtained on a 32^ x 8 lattice, together with the fit curves given by eq. ( |3.1| ). As one 
can see the position of the peaks decreases by increasing the external field strength. In 
Table 1 we report the values of the critical couplings versus next- 

Figure |5| displays instead the susceptibility of the absolute value of the Polyakov loop 
together with the susceptibility of the chiral condensate in the peak region for the largest 
explored value of the external magnetic field (next = 5). It is worthwhile to note that, as 
mentioned earlier, the position of the peak of Polyakov loop susceptibility (/3 = 5.2719(164)) 
and the position of the peak of the chiral condensate susceptibility (/? = 5.2694(84)) are 
consistent with each other and with the position of the peak obtained from the derivative 
of the free energy (/3 = 5.2680(38)), thus confirming the conclusion stated above, i.e. that 
the chiral and the deconfinement transition are shifted towards lower temperatures by the 
presence of the external field in an equal way, i.e. they continue to be coincident within 
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Figure 4: The derivative of the free energy eq. ( p.7[ ) versus (3 for some values of the strength of 
the constant chromomagnetic field parameterized (see eq. ( 2.11 )) by the integer ricxt- The yellow 
full circles corresponds to runs performed with different machines (APEmille) and algorithms as a 
check. 



lattice size 



next 



32 X 32 X 32 X 8 
64 X 32 X 32 X 8 
32 X 32 X 32 X 8 
32 X 32 X 32 X 8 
32 X 32 X 32 X 8 
32 X 32 X 32 X 8 
32 X 32 X 32 X 8 




1 
1 
2 
3 
4 
5 



5.4851 (202) 
5.4288 (128) 
5.3808 (128) 
5.3228 (90) 
5.2888 (44) 
5.2659 (48) 
5.2680 (38) 



Table 1: The values of the critical coupling versus the external field strengths. 



statistical errors even for next 7^ 0. 

The numerical results obtained so far let us conclude that the critical coupling is 
dependent on the strength of the background constant chromomagnetic field. On the other 
hand for pure SU(3) gauge theory we obtained [22] that the value of the critical coupling 
is not changed by a monopole background field. Similar results are expected in presence 
of dynamical fermions [16,23,24], however we will verify this fact explicitly for the present 
case. 

We recall the definition of an abelian monopole background field on the lattice, for 
more details and physical results see ref. [16]. It is well known that for SU(3) gauge theory 
the maximal abelian group is U(l)xU(l), therefore we may introduce two independent 
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Figure 5: The susceptibility of tlie absolute value of the Polyakov loop together with the suscep- 
tibility of the chiral condensate. The vertical full line represents the position of the peak in the 
derivative of the free energy for chromomagnetic field strength next — 5- The vertical dotted lines 
give the error region. Red and blue full lines are the best fits with the same parameterization as in 



eq. (3T) respectively to the susceptibility of the (absolute value of the) Polyakov loop and to that 
of the chiral condensate. 



types of abelian monopoles using respectively the Gell-Mann matrices A3 and Ag or their 
linear combinations. In the following we shall consider the abelian monopole field related 
to the A3 diagonal generator. In the continuum the abelian monopole field is given by 



a, 3 ^mon 



X X n 



\x\(\x\ 



X ■ n) 



(3.4) 



where n is the direction of the Dirac string and, according to the Dirac quantization 
condition, 

^moii QjH integer. The lattice links corresponding to the abelian monopole field 
eq. (^.41) are (we choose n = X3) 



^f|(f) 

ur{x) 



e 1.2 



e-'^'^r^^^ 
1 



(3.5) 



with e'^f'ix) defined as 



{X2 - X2) 



1 



+ - 



arir, 



l^^monl l^^n 

{xi - Xi) 



(X3 - X;) 
1 

(X3 - X3) 



(3.6) 
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where {Xi, X2, X^) are the monopole coordinates, Xmon = {x — X). The monopole back- 
ground field is introduced by constraining (see eq. ( |3.5| )) the spatial links exiting from the 
sites at the boundary of the time slice xt = 0. For what concern spatial links exiting from 
sites at the boundary of other time slices (xt 7^ 0) we constrain these links according to 
eq. Q. 

We have performed numerical simulations in presence of an abelian monopole back- 
ground field with monopole charge nmon = 10 (again for 2 staggered flavors QCD of mass 
aruq = 0.075). The critical coupling has been located by looking at the peak of the deriva- 
tive of the free energy (see figure (^)). We find 

= 5.4873(192) . (3.7) 



We have also done simulations in absence of any external chromomagnetic field, finding 
that the susceptibilities of the Polyakov loop, of the chiral condensate and of the plaquette 
display a peak at 

Pc = 5.495(25) . (3.8) 

Noticeably, this value of the critical coupling without external field is consistent, within 
our statistical uncertainty, with the value we get when we consider an abelian monopole 
field as background field. Therefore we can conclude that, as we found [16] in the case of 
pure gauge SU(3) gauge theory, the abelian monopole field has no effect on the position of 
the critical coupling. 
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4. Deconfinement temperature and critical field strength 



In previous studies [2] in pure lattice gauge theories we looked for the possible dependence 
of the deconfinement temperature on the strength of an external (chromo) magnetic field. 
In particular we studied SU(2), and SU(3) l.g.t.'s both in (2+1) and (3+1) dimensions and 
U(l) l.g.t. both in 4 dimensions and in (2+1) dimensions. In fact, in the case of non-abelian 
gauge theories, irrespective of the number of dimensions, we found that the deconfinement 
temperature depends on the strength of the constant chromomagnetic background field 
(similar studies have been performed within a different framework in refs. [25,26]). On the 
other hand, for U(l) gauge theory we found no evidence for a dependence of the critical 
coupling on the strength of an external magnetic field. In particular, as is well known, 4- 
dimcnsional U(l) l.g.t. undergoes a transition from a confined phase to a Coulomb phase: 
our analysis showed [2] that the location of the confinement- Coulomb phase transition is 
not changed by varying the strength of an applied constant magnetic field. The same 
analysis has been performed for compact quantum electrodynamics in (2+1) dimensions 
where it is known [27] that at zero temperature external charges are confined for all values 
of the coupling and it is well ascertained that the confining mechanism is the condensation 
of magnetic monopoles which gives rise to a linear confining potential and a non-zero string 
tension. Even in this case we verified that the critical temperature for deconfinement does 
not depend on the strength of an external magnetic field. As a consequence, we concluded 
that the dependence of the critical coupling on the strength of the external chromomagnetic 
field is a peculiar feature of non-abclian theories. 

The main aim of the present investigation is to extend our study to non-abelian gauge 
theories in presence of dynamical fcrmions. To this end, having determined in the previous 
section the critical couplings corresponding to different external field strengths, we now try 
to estimate the critical temperature: 

a[[ic-,mq)Lt 

where Lt is the lattice temporal size and a{f3c, ruq) is the lattice spacing at the given critical 
coupling (5c- In the case of SU(3) pure gauge theory in order to evaluate a(/3c) we used [2] 
the string tension. We obtained that the values of the critical temperature versus the 
square root of the external field strength can be fitted by a linear parameterization. By 
extrapolating to zero external field strength we obtained: 

= 0.643(15) (4.2) 

in very good agreement with the estimate Tc/yfa = 0.640(15) in the literature [28]. At 
the same time the intercept of line with the zero temperature axis furnished an estimate 
of the critical field strength (i.e. the limit value above which the gauge system is in the 
deconfined phase even at very low temperatures) 

^/gHl = (2.63 ± 0.15) = (1.104 ± 0.063) GeV = 6.26(2) x 10^^ Gauss (4.3) 
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using for the physical value of the string tension ^Ja = 420 MeV. The same analysis can 
be performed by means of the improved lattice scale introduced in Refs. [29,30] 



A = -f{g^){l + C2a{gf + c^a{gf) , a{g) 



fig' 



fig' 



(4.4) 



where g is the gauge coupling, C2 = 0.195(16), C4 = 0.0562(45), A/^/a = 0.0138(12) and 
fig"^) is the 2-loop scaling function 



fig') = ibog'r'^/''o 



exp 



1 



2bog' 



(4.5) 



with 



60 



1 



167r2 
1 



167r2 



Nc 2 
11— - -Nf 
3 3 ^ 



34 . /lO,, N. 
—N^ - —Nc + 
3 ' V 3 



N 



(4.6) 



Nc is the number of colors and Nf is the number of flavors. 

Obviously the analysis of pure gauge data in units of A gives results consistent with 
the same analysis done using the scale of the string tension (see Ref. [2]). In particular a 
linear extrapolation towards zero external field gives: 



A 



45.05(1.02) 



in agreement with Tc/A = 46.38(1.09) obtained from Tc 
critical field strength turns out to be 



(4.7) 



a = 0.640(15). Moreover the 



A 



209.6 ±3.07 



(4.8) 



which corresponds to ^JgHc = 1.21(ll)Gey in agreement with the estimate ^JgHc = 

1.10(6)Gey given in eq. (U). 

Let us turn now to the A'^^ = 2 case. Also in this case we have to face the problem of 

fixing the physical scale. In order to reduce the systematic effects involved in this procedure 

we will consider the ratios 

TcigH) _ ^ 



rj. VS. ^ (4.9) 

where Tc is the critical temperature without external field. The above quantities can be 
obtained once the ratio of the lattice spacings at the respective couplings is known. A 
rough estimate of this ratio can be inferred by using the 2-loop scaling function fig'^) given 



in eq. (|4.5| ) for Nf = 2. A better estimate could be obtained, as in the quenched case, 
by adopting an improved scaling function f(g'^){l + 020(5)^ + c^aig)^). We do not know, 
however, the values of C2 and C4 for Nf = 2. In a first approximation we will fix C2 and C4 
to their quenched values given above. In figure ^ the quantities reported in eq. ( [4. 91 ) are 



- 12 - 



^■25| I I I I I I I I I I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' 




VgH / 

Figure 7: The critical temperature TdgH) at a given strength of the chromomagnetic background 
field in units of the critical temperature Tc without external field versus the square root of the 
strength of the background field in the same units. Red circles are obtained by adopting the 
improved scaling function. The blue line is the linear best fit. The blue circle on the horizontal axis 
is the linear extrapolated values for the critical background field. Green squares are obtained by 
adopting the 2-loop scaling function. The blue dashed line is the linear best fit. The blue square 
on the horizontal axis is the linear extrapolated values for the critical field. 



displayed for both choices described above, i.e. 2-loop asymptotic scaling and improved 
scaling. 

The main result of our investigation is clear: even in presence of dynamical quarks the 
critical temperature decreases with the strength of the chromomagnetic field; moreover a 
linear fit to our data can be extrapolated to very low temperatures, leading to the prediction 
of a critical field strength above which strongly interacting matter should be deconfined at 
all temperatures. 

However, as can be appreciated from figure |^, the exact value of the critical field 
strength is largely dependent on the choice of the physical scale. In particular, assuming 
the 2-loop scaling function we obtain 

= 26.8(5) , (4.10) 
while assuming the improved scaling function we obtain 

^^^ = 4.29(10). (4.11) 

c 

If the deconfinement temperature at zero field strength is taken to be of the order of 170 
MeV, that means y/gHc in the range 0.7 — 4.5 GeV. 
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Figure 8: The chiral condensate eq. ( |3.2| ) versus /3 in correspondence of some values of the constant 
chromomagnetic background field. In the inset the region corresponding to the phase transition has 
been magnified. 



It is clear that in order to get a reliable estimate of the critical field strength, suitable 
for phenomenological purposes, one should get a more reliable estimate of the physical 
scale of our lattices. However one should work at a fixed value of the pion mass (as 
close as possible to its physical value) as well: that is beyond the purpose of the present 
investigation and will be the subject of further studies. 

To conclude this section we consider our measurements of the chiral condensate. In 
figure |8| we display the chiral condensate versus the gauge coupling in correspondence of 
some values of the external field strength. Our numerical data show that, at least in the 
critical region, the value of the chiral condensate depends on the strength of the applied 
field. Interestingly enough, similar results for the chiral condensate have been found in 
ref. [31]. 



5. Summary and Conclusions 

In this paper we have studied how a constant chromomagnetic field perturbs the QCD 
dynamics. In particular we focused on the theory at finite temperature and we have found 
that, analogously to what happens in the pure gauge theory [2], the critical temperature 
depends on the strength of the constant chromomagnetic background field: it decreases as 
the external field is increased and we have inferred, as an extrapolation of our results, that 
eventually the system is always deconfined for strong enough field strengths. We estimated 
this critical field strength to be of the order of 1 GeV, which is a typical QCD scale [32]. 
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Notice that our estimate is of course affected by several systematic uncertainties, like that 
regarding the estimate of the physical scale: for that reason we consider it only as an order 
of magnitude of the real expected critical strength. In fact a more reliable determination 
usable for phenomenological purposes should be performed by working with a fixed physical 
value of the pion mass (and as close as possible to its physical value); that is out of the 
purpose of the present work and will be the subject of future investigations. 

By comparing the critical couplings determined from the derivative of the free energy 
functional with those determined from the susceptibility of the chiral condensate and of 
the Polyakov loop we have ascertained that, even in presence of an external chromomag- 
netic background field and at least up to the field strengths explored in the present work, 
the critical temperatures where deconfinement and chiral symmetry restoration take place 
coincide within errors. 

Another intriguing aspect we have found is the dependence of the chiral condensate on 
the chromomagnetic field strength. This last point deserves further studies. In order to get 
a deeper understanding of our results, we also plan to study the effect of the background 
field on the equation of state of QCD. 
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